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1 .0  INTRODUCTION 
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In  previous  reports  on  this  topic  *  ,  we  have  reviewed  progress  in  our 
study  of  the  three-dimensional  boundary  layers  on  bodies  of  revolution  when 
the  external  velocity  is  prescribed.  The  governing  equations  are  parabolic  in 
y,  the  coordinate  normal  to  the  body,  and  x,  the  direction  of  the  local  veloc¬ 
ity  parallel  to  the  body,  and  so,  were  it  not  for  two  features  of  the  problem, 
standard  methods,  which  have  proved  to  be  successful  in  two  dimensions,  might 
readily  be  adapted  to  integrate  them.  The  first  feature,  which  largely  formed 
the  topic  of  the  first  report^,  is  that  the  obvious  body-based  coordinate 
system  leads  to  equations  which  contain  artificial  singularities.  An  approp¬ 
riate  coordinate-transformation  is  needed  to  remove  them  and  this  requires 
some  care  to  devise  but  once  it  is  done  the  coefficients  of  the  various  deriv¬ 
atives  in  the  equations  are  bounded  everywhere  and  they  are  well -adapted  to 
standard  numerical  procedures. 

The  second  feature  is  that  the  local  streamwise  direction  in  the  boundary 
layer  varies  from  point  to  point  both  laterally  and  longitudinally.  Eventu¬ 
ally  as  separation  is  approached,  especially  on  the  leeside,  this  gives  rise 
to  difficulties  as  it  may  be  unavoidable  that,  over  part  of  the  boundary- 
layer,  Integration  proceeds  against  the  local  flow  direction..  Thus  it  is 
usual  to  compute  the  solution  on  the  windward  (<t>  *  0)  and  leeward  (<j>  »  tt) 
line  of  symmetries  first  and  then  at  each  x-station  (where  x  measures  distance 
along  the  body  from  the  nose)  Integrate  over  the  body  either  from  <J>  *  0  to 
<t>  *  it  or  vice-versa.  For  such  a  procedure  to  be  successful  using  standard 
methods,  it  is  Important  that  w  (the  velocity  component  in  the  4>  direction) 
should  not  change  sign  as  y  varies  for  fixed  x,<j>.  If  it  does,  then  which¬ 
ever  direction  of  <)>  we  Integrate  in,  over  part  of  the  range  of  y  we  must  be 
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Integrating  against  the  local  stream  direction  and  hence  not  making  use  of  all 

2 

the  Information  needed  to  compute  the  solution.  In  the  second  report  a 

remedy  for  this  problem  was  devised.  The  numerical  method  in  straightforward 

computation  of  boundary  layers  used  by  us  is  the  standard  box  originally 
3 

devised  by  Keller  .  Once  w  changes  sign,  we  replace  it  either  by  the  zigzag 
4  5 

box  or  by  the  characteristic  box'.  Both  of  these  methods  are  an  Improve¬ 
ment  on  the  standard  box  In  such  circumstances  because  they  take  account  of 

4 

the  extra  physical  Information  available.  The  zig-zag  box  Is  simpler  to 
use  but  is  less  sensitive  and  produces  some  undesirable  features  In  the 
computed  solution  as  the  limit  of  computability  Is  reached.  The  charac- 

5 

terlstic  box  is  more  complicated  to  program  but  seems  largely  free  of  spur¬ 
ious  features  and  is  able  to  compute  the  solution  practically  up  to  the  mathe¬ 
matical  limits  of  the  formulation. 

This  limit  Is  the  accessibility  boundary  and  there  are  precise  mathemati¬ 
cal  reasons  why  the  computation  cannot  be  extended  beyond  It.  Partly  it  con¬ 
sists  of  the  separation  line,  at  which  the  solution  develops  a  singularity  and 
partly  by  an  external  streamline  passing  through  the  forward  separation  point 
A,  known  as  the  ok  (Turk,  arrow)  of  accessibility.  The  situation  is  Illus¬ 
trated  In  Fig.  1. 


*  =  it  LEEWARD  LINE  OF  SYMMETRY 


...  *  =  0  WINDWARD  LINE  OF  SYMMETRY 

Fig.  1.  The  solution  can  be  computed  over  the  hatched  part  of  the  body  start¬ 
ing  from  the  forward  stagnation  point.  Notice  that  the  nose  and 
tails  are  single  points  but  In  polar  coordinates  they  must  be 
replaced  by  lines  (0  <  <p  <  it). 
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Future  work  envisaged  at  the  completion  of  these  studies  was  to  interact 
the  boundary-layer  computation  with  an  inviscid  code  which  would  thereupon 
remove  the  barrier  imposed  by  the  singularity  at  separation  and  enable  us  to 
compute  the  flow  field  over  the  whole  body.  We  also  wished  to  study  other 
body  shapes  to  investigate  how  far  the  conclusions  from  this  work  could  be 
regarded  as  general. 

The  first  part  of  such  work  is  a  large  program  of  study  and  has  proceeded 
independently  of  the  work  to  be  described  in  this  report  but  it  is  useful  to 
comment  briefly  on  what  has  been  achieved.  Interactive  studies  have  been  com¬ 
pleted  on  two  problems  in  two-dimensional  flow6,7.  We  have  examined  the 
boundary  layer  near  the  leading  edge  of  a  thin  airfoil  at  incidence6,7  and 
shown  that  once  interaction  with  the  external  stream  is  permitted,  the  exist¬ 
ence  of  separation  no  longer  prevents  the  successful  completion  of  the  solu- 

o 

tion.  Comparison  with  experimental  data  obtained  by  Gault  is  generally 

g 

favorable  .  However,  there  is  a  curious  feature  of  the  solution  which 
appears  to  be  a  genuine  property  of  high-Reynolds  number  flows.  This  is  that 
although  reversed  flow  can  be  found  but  its  extent  is  limited  before  the  whole 
iteration  procedure  fails  to  converge.  A  mathematical  analysis  of  the 
solution7  confirms  the  correctness  of  the  general  conclusion  and  it  is  also 
broadly  unaltered  if  transition  to  a  turbulent  boundary  layer  occurs,  the 
turbulence  being  modelled  by  an  algebraic  eddy  viscosity.  It  appears  that  the 
notion  of  a  thin  boundary  layer  on  the  airfoil  modifying  an  external  potential 
flow  has  a  certain  range  of  validity  in  which  is  included  the  occurrence  of  a 
small  region  of  reversed  flow  in  the  boundary  layer,  but  thereafter  becomes 
Inappropriate.  One  may  crudely  regard  the  approach  as  going  sour.  In  view  of 
this,  it  was  decided  to  carry  out  a  further  comparison,  with  Briley's  solu¬ 
tion,  using  the  Navier-Stokes  equations  of  a  problem  with  separation10. 


3 


Again  good  agreement  was  obtained  with  his  published  data  and  again  the 
boundary-layer  approach  eventually  went  sour.  In  the  same  conditions, 

Briley's  numerical  procedure  failed  to  converge.  Although  he  attributed  this 
to  a  physical  instability  of  the  flow,  we  have  reservations  and  believe  that 
dramatic  changes  are  occurring  In  the  solution  for  which  his  domain  of  inte¬ 
gration  is  too  small  to  describe  adequately.  Work  is  in  progress  to  explore 
this  problem  further. 

The  subject  of  this  report  is  the  progress  made  in  applying  boundary-layer 
theory,  without  interaction,  to  the  computation  of  supersonic  flow  past  cones 
at  incidence.  In  its  simplest  form  this  is  a  two-dimensional  problem  because 
the  streamwise  coordinate  can  be  absorbed  in  the  formulation  and  the  independ¬ 
ent  variables  are  n,  the  reduced  normal  distance  from  the  cone  and  <t>,  the 
azimuthal  angle  measured  from  the  windward  line  of  symmetry.  Solutions  can 
always  be  found  at  <J>  *  0  and  continued  round  towards  <j>  =  tt,  the  leeward 
line  of  symmetry,  but  are  difficult  to  find  near  4>  =  u  even  if  separation 
does  not  occur  beforehand.  At  <|>  *  it  pairs  of  solutions  can  be  found  in 
certain  circumstances  but  none  in  others  and  the  situation  is  not  always  com¬ 
patible  with  the  solution  properties  for  <|>  <  tt.  One  aim  of  this  report  is 
to  reconcile  these  features  of  the  flow. 

Commonly,  cones  do  not  occur  alone  in  practical  configurations  but  are 
part  of  more  complicated  shapes  having  attached  to  them  fins  and  afterbodies, 
which  may  be  boattailed  and  asymmetric,  and  possibly  rounded  noses.  It  was 
thought  appropriate,  therefore,  to  set  the  semi-similar  solutions,  described 
above,  in  a  more  general  context  and  with  this  in  mind  we  introduced  suction 
near  the  vertex  of  the  cone  so  that  a  self-similar  solution  is  only  generated 
at  some  distance  along  it.  We  expect  that  if  separation  occurs,  there  would 
also  exist  an  accessibility  boundary  and  we  should  be  able  to  Investigate  its 
properties.  Looking  further  ahead,  we  shall  then  be  able  to  comment  on  the 


parabolic  elliptic  boundary  layer  (PEBL)  method  of  computing  the  boundary 

layer  once  separation  has  occurred.  We  have  reservations  about  this  method 

since  there  is  a  nonphysical  basis  to  it.  Another,  more  promising,  method  is 

11  12 

the  use  of  parabolized  Navier-Stokes  equations  *  and  we  expect  that  the 
extension  of  our  work  to  include  interaction  will  lead  to  a  useful  alternative 
way  of  tackling  practical  problems  and  have  a  status  similar  to  that  in  the 
two-dimensional  flows  discussed  earlier. 


2.0  GOVERNING  EQUATIONS 


The  governing  boundary-layer  equations  for  compressible  laminar  and  turbu¬ 
lent  flows  on  a  general  conical  surface  of  semi-vertex  angle  0C  are  well 
known  and  can  be  written  in  the  following  form: 

Continuity: 


lx  (pux)  +  te  (pw)  +  |y  (pVX>  =  0 


(2.1) 


x-Momentum 
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Energy 


,U1H  +ew  SH  t-3H 
3X  X  36  3y 


3  fy  3H 
3y  |_Pr  3y 


+  u(1  -fc>  k  -PvnrQ(2.4) 


Here  p7  =  pv  +  p ' v ' ;  0  denotes  the  polar  coordinate  in  the  developed 
plane;  x  the  coordinate  along  the  generators;  and  y  the  coordinate  normal  to 
the  surface;  with  w,  u,  v  the  velocities  in  the  0,x  and  y  directions. 

At  the  windward  and  leeward  stagnation  lines  where  0=0  and  tt  sinBc,  the 
cross-flow  momentum  equation  is  an  identity  since  w  =  dp/d0  =  0.  Taking 
into  account  the  symmetry  conditions  and  differentiating  Eq.  (2.3)  with 
respect  to  0,  we  can  write  it  as 
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where  wQ  =  3w/30.  Noting  that  w  =  0,  Eqs.  (2.1),  (2.2)  and  (2.4)  can 
be  simplified  and  written  as 

fj  (pux)  +  Pw0  +  |p-  (pvx)  =  0 


(2.6) 
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pu  fx  +  pv  ly  =  |y  [fcly  +  u  0  “Jr*  ly  "  pv  H  ] 


(2.7) 


(2.8) 


Before  we  discuss  the  solution  of  the  above  equations,  we  express  them  in 
terms  of  similarity  variables.  For  the  general  case  in  which  the  governing 
equations  are  given  by  those  defined  by  Eqs.  (2.1)  to  (2.4),  we  define 


u. 


1/2 


dr'-(vv)  ody 


In  addition,  we  Introduce  a  two-component  vector  potential  such  that 


(2.9) 
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l^x1  /2g(x,e,n) 

(2.11b) 

Using  these  transformations  and  the  eddy  viscosity  and  turbulent  Prandtl 
number  concepts,  that  Is, 
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We  can  write  Eqs.  (2.1)  to  (2.4)  as 
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Here  primes  denote  differentiation  with  respect  to  n  and 

b  .  c(i  +  e;i 

(2-16) 
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Similarly  the  line-of -symmetry  equations,  Eqs.  (2.5)  to  (2.8)  can  be  trans¬ 
formed  and  expressed  in  a  form  similar  to  those  given  by  Eqs.  (2.13)  to  (2.15) 
by  using  the  transformations  given  by  Eqs.  (2.9)  to  (2.11)  with  slight  modifi¬ 
cation  to  the  following  parameters 
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Equations  (2.5),  (2.7)  and  (2.8)  then  become 
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b,E‘  +  j,-®-  0  F'F'j  ♦<§F+a  +  J>)E'  =  x(F'  |f -  £•  |f)  (2-2') 

i  e 

i 

Here  6'  3  wQ/Ug. 

The  boundary  conditions  are 

i 

f‘  »  g‘  *  g  ■  F'  -  B'  •  6  *  0  at  r,  =  0  (2.22) 


and 


f '  =  E  =  F '  ■  1, 


g'  =  G' 


ua  d0 
e 


at  n  =  ro 


In  addition  either  E'  or  E  may  be  prescribed  at  n  =  0  and  in  order  to  start  the 
computations  we  shall  apply  suction  to  the  boundary  layer  in  0  <  x  <  0.05. 

Thus  f  is  a  prescribed  function  of  0  when  n  *  0  for  these  values  of  x 
and  thereafter  f  *  0  at  the  wall.  It  is  convenient  to  introduce  a  parameter  k 
in  the  line-of -symmetry  studies  defined  by  Moore  to  be 


k  3 

The  explicit  forms  for  ug 
theory  and  are 


,2d,e) 

and  we  follow  from  inviscid  conical  flow 


(2.23) 


ue  3  1  +  A  sinec  -  A  s1n&c  cos(e/sinec) ,  wg  3  A  sin(e/sinec)  (2.24) 

where  0C  is  the  semi -vertex  angle  of  the  cone  and  A  is  a  constant  depending  on 
Bc,  and  a,  the  angle  of  attack  of  the  cone.  Typically  A  %  1/4  and 
explicit  values  have  been  tabulated  by  Sims14.  The  numerical  calculations 
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to  be  described  in  Sec.  3  were  carried  out  for  £c  *  10  .  It  will  be 
convenient  subsequently  to  define  4>  *  9/sin3e,  to  be  the  angular 
distance  around  the  cone  from  the  windward  line  of  symmetry  so  that  <t>  =  tt 
is  the  leeward  line  of  symmetry. 

Since  the  studies  are  basic  in  nature,  there  is  no  advantage  in  solving 


the  governing  equations  for  a  general  fluid.  Consequently  we  set 


Pr  *  1  and  up  ■  u^>e»  i.e.  C=1  (2.25) 

A  further  simplification  is  to  assume  that  the  thermal  boundary  condition  is 
either  that  cone  is  adiabatic  or  the  wall  temperature  is  prescribed.  Then 
(2.15)  reduces  to 


i.e.  C  =  1 


(2.25] 


E  =  Ew(l  -  S)  +  S 


(2.26) 


where  S  =  1  or  F',  respectively,  for  the  two  cases.  Hence 


-  0  +|’(y-1)M^)(S  +  Ew(1-S))-]-(y-1)^F'2  (2-27) 


3.0  RESULTS 


We  have  used  the  box  method  to  solve  the  governing  equations  discussed  in 
the  previous  section.  In  regions  where  the  circumferential  velocity  component 
changed  sign  becoming  negative,  we  used  the  zig-zag  box  method.  Since  both  of 
these  methods  are  already  discussed  in  detail  in  several  references^’  , 
we  shall  not  elaborate  them  here  but  only  present  the  results  of  the 
computations. 

3.1  The  Leeside  Line  of  Symmetry  in  Laminar  Flow 

Since  we  have  an  extensive  body  of  new  results  for  this  aspect  of  the  flow 
field,  it  is  worthwhile  to  begin  our  discussion  of  them  with  a  review  of  pre¬ 
vious  work. 

The  boundary-layer  equations  appropriate  to  the  windward  and  leeward  lines 
of  symmetry  of  a  cone  fixed  in  a  supersonic  external  stream  were  first  written 
down  by  Moore  in  terms  of  a  parameter  k  (positive  for  the  windward  line 
of  symmetry  and  negative  for  the  leeward,  defined  in  (2.23)).  The  value  of 
k,  for  a  prescribed  external  Mach  number  and  ratio  of  wall  to  freestream 
enthalpy,  determines  the  existence  or  otherwise  of  a  solution  of  these  simi¬ 
larity  equations.  The  solution  of  these  equations  in  our  independent  variable 
exist  for  all  k  <  0  and  for  certain  k  >  0  though  they  may  not  be  unique. 
They  are  expected  to  be  solutions  of  the  full  boundary-layer  equations  in  the 
limit  x  ♦  “  where  x  measures  distance  along  a  generator,  and  as  4>  0  or  n 

where  $  -  0  is  the  windward  plane  of  symmetry  and  <j>  =  n  the  leeward.  These 

16  17  18  19 

equations  have  been  considered  by  Cheng  ,  Roux  ,  Murdock  ,  Wu  and  Libby 

20 

and  Rubin,  Line  and  Tarulli  . 


One  way  of  attempting  to  gain  an  understanding  of  the  nonuniqueness  and 
nonexistence  of  these  similarity  solutions  is  to  explore  the  independent  vari¬ 
able  x  and  to  integrate  the  nonsimilar  equations  for  the  symmetry  planes 

along  the  generator  with  increasing  distance  from  the  vertex  as  was  undertaken 

20 

for  example  by  Rubin  et  al  .  The  solution  so  found  either  breaks  down  at  a 
finite  value  of  x,  or  it  may  be  continued  to  x  «  »  the  limit  solution  pos¬ 
sibly  being  given  by  the  similarity  equations.  Another  way  is  to  restore  <|> 

rather  than  x  and  to  perform  a  windward  to  leeward  integration  in  the  manner 
21  22  17  23 

of  Cooke-  ,  Boericke  ,  Roux  ,  and  Tarulli  .  If  <j>  =  it  is  attained  then  the 
limit  solution  may  or  may  not  be  a  solution  of  the  similarity  equations.  A 

22 

third  approach  is  to  solve  the  complete  equations  in  the  manner  of  Boericke 
24 

and  Dwyer  .  Since  there  are  three  parameters  involved  and  the  similarity 

equations  have  nonunique  solutions,  a  complete  documentation  of  the  relation 

behavior  the  various  approaches  has  not  yet  been  achieved. 

One  of  the  aims  of  the  present  investigation  is  to  examine  the  possibility 

that  the  solution  of  the  nonsimilar  equations  for  the  leeward  plane  of  sym- 

20 

metry  comes  to  an  end  at  a  finite  value  of  x.  Rubin  et  al.  noted  that 
for  those  k  for  which  they  encountered  this  phenomenon,  the  streamwise  skin 
friction  vanished  though  this  did  not  seem  to  be  a  separation  of  the  Goldstein 
type  in  that  it  was  caused  not  by  an  adverse  pressure  gradient  but  by  a  viola¬ 
tion  of  the  assumption  that  the  crossflow  is  zero  on  the  plane  of  symmetry. 

The  crossflow  boundary  layers  were  colliding  and  forcing  the  streamwise  bound¬ 
ary  layer  to  leave  the  wall.  The  situation  is  very  similar  to  that  encount- 

?5 

ered  by  Stewartson,  Cebeci  and  Chang  in  their  study  of  the  development  of 
a  laminar  boundary  layer  near  the  entrance  of  a  loosely  coiled  pipe.  On  the 
inner  line  of  symmetry  there  appeared  to  be  a  local  singularity  where  the 
axial  component  of  skin  friction  vanished.  Further  downstream,  the  flow  was 


well-behaved  except  in  the  Immediate  neighborhood  of  the  inner  generator.  The 

two  secondary  boundary  layers  had  collided.  The  behavior  of  the  solution  near 

the  local  singularity  has  more  recently  been  analyzed  by  Stewartson  and 
26 

Simpson  and  our  aim  in  this  study  is  to  demonstrate  that  their  analysis  is 
also  applicable  to  the  colliding  boundary  layers  on  the  leeward  side  of  the 
cone. 

We  shall  show  that  a  singularity  is  theoretically  possible  if  0  >  k  >  -2/3, 

regardless  of  the  Mach  number  Mg  and  enthalpy  ratio  Ew,  and,  further,  that 

much  of  the  analysis  of  ref.  26  may  be  taken  over  without  revision  although 

there  are  significant  differences  in  the  outer  part  of  the  terminal  boundary 

layer.  Similarity  solutions  are  possible  in  two  ranges  of  values  of  k  and 

then,  most  probably,  occur  in  pairs.  The  first  range  is  k  >  k^*  where 

k^  >  -2/3  is  a  number  depending  on  the  Mach  number  Mg  and  the  enthalpy  ratio 

Ew.  Thus,  if  Mg  =  0  we  have  found  numerically  that  k-j*  *  -0.292  but  no 

attempt  was  made  to  find  any  dual  solutions  in  k  >  k-|*.  The  existence  of  the 

1  ft 

second  solution  was  established  by  Murdock  In  computations  with  Ew  =  0.4,  and 

p 

Mg  *  1 6/(y  -  1),  when  k^  *  -0.079.  A  second  range  of  k  for  which 
similarity  solutions  are  possible  is  k2*  >  k  >  -1  where  k?*  is  very  close 
to  -2/3  but  distinct  from  it.  Thus  Murdock  gives  k2*  a  -0.6656.  If  Mg  ®  0  we 
were  unable  to  find  a  similarity  solution  when  k  *  -0.66  and  when  k  «  -2/3 
there  are  at  least  two  solutions  for  all  Mg,  E^3,1®. 

Numerical  solutions  of  the  nonsimilarity  leeside  equations  are  found  here 
for  a  number  of  values  of  k  and  Mfi  *  0,  it  being  considered  that  the  charac¬ 
ter  of  the  solutions  Is  Independent  of  Mg,  Ew.  If  k  *  -0.25  the  solution 
may  be  found  for  all  distances  x  along  the  leeside  generator  and  it  approaches 
the  primary  similarity  solution  as  x  +  «,  and  this  is  also  the  limit  of  the 


corresponding  equations  In  which  the  x-varlatlon  Is  suppressed  and  <j>  allowed 
to  Increase  from  $  *  0  at  the  windward  line  of  symmetry  to  $  *  n.  If  k  *  -0.5 
the  solution  terminated  at  a  finite  value  of  x  In  a  singularity  In  accord  with 
our  theoretical  considerations.  If  k  *  -0.70  It  appears  that  the  solution  may 
be  continued  Indefinitely  in  x.  We  Infer  that  a  singularity  will  occur  If 
k-j*  >  k  >  k2*. 

3.2  Numerical  Studies 

Equations  (2.19)  to  (2.21)  hold  for  either  the  windward  or  leeward  planes 
of  symmetry.  At  the  windward  plane  of  symmetry  k  >  0  and  at  the  leeward 
plane  k  <  0.  There  are  three  nondlmenslonal  parameters  in  the  problem, 
namely  Mg,  k  and  Ew,  and  both  the  nonslmilar  and  similar  forms  of  (2.19)  to 

(2.21)  have  been  considered  by  various  authors.  A  useful  summary  of  earlier 

20 

work  is  given  by  Rubin,  Lin  and  Tarulli  .  Dual  similarity  solutions  exist 

for  k  >  k.|*  and  again  when  -1  <  k  <  k2*  where  k2*  &  k^*  <  0  and  k^*,  k2* 

depend  on  the  values  of  Mg  and  Ew  (see  refs.  16  to  19).  When  k^*  >  k  >  k2*  no 

21 

solutions  of  the  similarity  equations  exist  and  solutions  by  Cooke  ,  Boericke 

among  others,  from  the  windward  plane  of  symmetry  to  the  leeward  plane  of 

symmetry  of  the  equations  for  the  cone  with  the  4>  dependence  retained  but 

with  x  ■  «>  indicate  that,  for  such  k,  the  assumption  implicit  in  (2.19)  to 

(2.21)  that  the  azimuthal  velocity  Is  proportional  to  n  -  <j>  Is  Invalid. 

Rather  it  remains  nonzero  as  $  tt  with  the  result  that  the  boundary  layers 

growing  on  either  side  of  the  cone  will  collide.  Integration  of  (2.19), 

20 

(2.20)  for  certain  k  In  this  range  by  Rubin  et  al.  demonstrates  the 
existence  of  a  point  xg  at  which  the  streamwlse  component  of  skin  friction 
(Uy  at  y  *  0)  vanishes.  As  the  Incidence  Is  increased  so  that  k  is  less 
than  k2*  It  Is  likely  that  there  is  cross-flow  reversal  In  the  windward  to 
leeward  Integration  before  <j>  reaches  tt. 


It  is  the  collision  phenomenon  and  the  existence  of  the  "separation"  point 

on 

xg  of  Rubin  et  al.  that  interests  us  here.  As  this  appears  regardless 
of  the  values  of  Ma  and  E^  we  set  M  =  0  and  Ew  s  1  in  (2.19)  to  (2.21)  and 
found  solutions  of  these  equations  and  of  the  equations  considered  by 
Boericke  and  Dwyer*  for  negative  values  of  k.  The  integration  of  the 
leeside  equations  is  initiated  by  assuming  that  the  cone  is  permeable  in  0  <_  x 
<  0.05  and  that  F(P)  >  0  in  this  range,  but  F(P)  -*■  0  as  x  -*■  0.05.  The 
precise  form  of  P  is  given  in  (3.26)  below.  The  initial  profile  is  then 
computed  at  x  =  0,  evolves  as  x  increases,  and  either  terminates  at  a  finite 
value  of  x  in  a  collision  or  may  be  continued  to  x  =  ®.  The  only  noteworthy 
characteristic  of  the  solution  is  that  the  boundary-layer  thickness  increases 
rapidly  near  any  breakdown  point  but  it  was  then  only  found  necessary  to 
increase  the  number  of  points  across  the  boundary  layer.  The  step  length  in  n 
is  0.1  and  that  in  x  is  in  the  range  0.0015  to  0.025,  consistent  with  accuracy. 
The  numerical  method  was  also  used  to  solve  the  similarity  equations  ((2.19) 
with  3/3x  *  0)  and  also  the  windward-leeward  equations  the  step  length  in 
<t>,  where  <J>  sin3c  *  0,  varying  from  5°  near  <j>  s  0  to  as  little  as 
0.1°  near  <j>  *  180°. 

The  results  are  in  general  agreement  with  those  of  previous  workers.  For 
the  similarity  equation  we  find  that  kj*  *  -0.292  and  k?*  *  -0.666.  In 
Fig.  2a  we  display  the  wall  shear  parameter  F"w  computed  from  the  similarity 
equations  and  from  the  <J>-wise  Integration.  In  the  range  kj*  <  k  <  0  the 
primary  similarity  solution  is  also  the  limit  solution  of  the  ^-integration.  For 
|  k^*  >  k  >  kg*  where  k3*  »  -0.80,  the  <j>- integration  leads  to  a  collision  at 

<J)  *  180°  as  can  also  be  seen  from  Fig.  2b  in  which  we  plot  the  value  of  (it  - 

II 

4>)G^  In  the  limit  <j>  *  it.  The  similarity  solutions  for  k  <  k2*  do 
1  not  appear  to  be  relevant  to  the  (^-integration  solutions.  For  k  <  k3* 
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separation  occurs  before  the  leeside  is  reached.  In  Fig.  3  we  display  the 
results  of  the  leeside  integrations.  For  k  =  -0.25  the  solution  appears  to 
exist  for  all  x  and  to  approach  the  primary  similarity  solution  as  x  +  «. 

II 

For  k  =  -0.50  the  solution  exhibits  the  "separation"  point  x$  with  F  . 

For  k  *  -0.65,  -0.66  no  initial  solutions  could  be  found  for  any  values  of 
P0  =  3/2  F (0,0) .  For  k  =  -0.70  the  solution  can  apparently  be  continued  for 
all  x,  and  approaches  one  of  the  similarity  solutions  as  x  »  is  PQ  < 

0.40.  The  fate  of  the  solution  with  PQ  =  0.70  as  x  -*>  =>  is  not  clear. 

We  shall  not  adopt  the  term  separation  for  the  point  x$  as  it  is  not  the 
usual  Goldstein  form  but  is  reminiscent  of  the  heralding  of  the  collision  of 
the  boundary  layers  on  the  line  of  symmetry  of  a  curved  pipe  recently  studied 

pc  pc 

by  Stewartson,  Cebeci  and  Chang  and  by  Stewartson  and  Simpson  .  In  the 
latter  paper  the  authors  showed  that  the  streamwise  skin  friction  vanishes 
with  a  singularity  that  is  not  of  the  Goldstein  type,  but  is  weaker,  and  sig¬ 
nals  a  lifting  of  the  streamwise  boundary  layer  from  the  wall  as  the  crossflow 
boundary  layers  presumably  collide  underneath  it.  They  present  an  analysis 
which  agrees  well  with  their  numerical  study.  We  now  show  that  if 
0  >  k  >  -2/3,  the  same  analysis  may  be  applied  to  equations  (2.19)  - 
(2.21)  with  encouraging  agreement  with  a  representative  computation  carried 
out  when  k  =  -0.5  and  M  *  0,  E  *  1,  the  incompressible  case. 

C  W 

3.3  The  Singularity  at  the  Start  of  Collision  on  the  Leeward  Side  of  Symmetry 
The  theory  of  ref.  26  for  the  curved  pipe  may  be  applied  directly  to  the 
system  (2.19)  -  (2.21)  for  all  Mg  and  Ew  provided  3/2  k2  +  k  >  0.  In 
ref.  26  the  authors  assume  the  existence  of  a  point  x  ■  x$  at  which  the 
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0.6 


Fig.  2a.  The  wall  shear  parameter  FMW  from  the  similarity  equations  (contin¬ 
uous  curve)  and  the  $-wise  Integration  (crosses).  S  denotes  the  onset 
of  separation  of  the  crossflow  before  $  reached  ir. 


Fig.  2b.  The  values  of  10*  (it  -  $)GMy  as  $  -*•  it  In  the  $-w1se  Inte¬ 
gration.  S  denotes  the  onset  of" separation  of  the  crossflow  for 
♦  <  ir  and  E  denotes  k  «  kt>  where  the  similarity  solution  comes  to 
and  end.  1 


k  =  - 


=  :8:?B 


k  =  -0.25 


k  =  -0.70 


A 

Fig.  3b.  The _ values  of  G"# from  the  leeslde  Integration.  The  arrows  Indi¬ 
cate  the  similarity  limits  for  k  =  -0.25  and  k  =  -0.70.  The 
plot  starts  at  x  «  0.05  where  the  suction  was  terminated. 
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streamwise  skin  friction,  F"  ,,  vanishes  and  show  that  a  consistent  structure 

w 

may  be  obtained  in  its  neighborhood.  In  order  to  facilitate  comparison  with 
their  work  we  first  define 

x  =  eX,  o  =  K_1/4Y,  F'  =  K1/2U,  G'  =  K1/2  W,  -  |  F  +  G  -  x  =  K1/4V 

(3.1) 

where  K  is  a  positive  constant  which  later  on  we  shall  choose  to  be 

-|k  (f  k  +  i)pe/pw 


where  pw  is  the  value  of  p  on  the  wall.  Then  equations  (2.19)  may  be 
replaced  by 


„  3U  .  »  3U  _  a2U 

U  ax  +  V 


u  +  V  —  +  UW  -W2+— —  k  +  1 )  —  = 

u  ax  V  aY  uw  w  2  K  l2  K  p  Y2 

iR+iX  +  1  u  .  K  «  o 
ax  3Y  2  u  w  u 


(3.2a) 


(3.2b) 

(3.2c) 


Also  (2.22),  (2.23)  become 

U  =  V  =  W  =  0  on  Y  =  0 


(3.3) 


U  -*•  K"1/2,  W -3/2  kK"1/2  as  Y  -►  it  (3.4) 

together  with 

■“»  [1  +  \  (Y  -  1)M2][S  +  Ew(l  -  S)3  -  1  (y  -  1)M2KU2  (3.5) 

where  S  =  1  if  the  wall  is  adiabatic,  and  S  =  K^2U  if  the  wall  temperature 
is  prescribed.  In  the  above,  U  is  simply  a  scaled  velocity  component  along 
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the  leeward  line  of  symmetry,  W  a  scaled  azimuthal  velocity  component  with  the 
factor  it  -  4>  removed,  and  V  is  related  to  the  normal  velocity. 


In  ref.  26  the  authors  showed  that  the  flow  in  the  neighborhood  of  the 
singularity  required  analysis  in  three  different  regions.  In  the  inner 
region  U  a  (X$  ■  X)  and  is  negligible  in  (3.2),  (3.5)  except  where  the 
term  also  contains  the  derivative  with  respect  to  X  as  these  are  assumed  to  be 
large.  Thus  to  leading  order  in  Xg  -  X  we  may  ignore  the  term  UW  in  (3.2b), 
the  term  3/2  U  in  (3.2c)  and  (3.5)  becomes  p  =  p  .  We  now  choose 
K  =  -3/2  k  (3/2  k  +  1 )p  /p  and  then  in  the  inner  layer  the  equations 

c  W 

are  exactly  those  that  hold  in  the  neighborhood  of  the  inner  generator  of  the 
curved  pipe  of  ref.  26  when  terms  of  relative  order  X$  -  X  are  neglected. 

Then,  as  there,  the  solution  in  the  inner  layer  is  of  the  form 

CO 

U  =  y2/3(X$  -  X){n-|  +  £  p”4n/3un(ni )>  +  0{(XS  -  X)2}, 

n=l 

co 

V  =  y]/3{ni  +  £  y-4n/3vn(ni)}  +  o(Xs  -  X)  (3,r>l 

n=l 


W  -  J  y“4n/3Wn(ni )}  +  0(XS  -  X) 

n=l 


the  only  difference,  apart  from  the  different  usage  for  the  symbols  U,  V,  W, 
(U,  V,  W  of  ref.  26  correspond  to  V,  W,  U  here),  being  the  order  of  the  error 
terms  in  (3.6).  In  the  above  y^»l  and 


nl 


'  ala2 


4 /3 

log(p1/  +  a-jo^)  =  —  log 


(3.7) 


where  *  0.5325,  a2  =  -0.4636  are  constants  determined  by  the  singularity 
itself  and  B-j  is  a  constant  that  depends  on  the  previous  history  of  the 


flow  for  X  <  Xs.  In  ref.  26  terms  of  (3.6)  up  to  n  *  2  were  explicitly 
calculated,  and  they  are  exactly  the  same  in  this  case.  All  we  shall  require 
here  in  detail  is  the  behavior  of  U  and  W  of  (3.6)  when  n-j  »  1  as  this 
is  needed  for  a  match  with  the  outer  solution.  These  are  obtained  directly 
from  (4.11)  of  ref.  26  as 


U  %  (Xs  -  X)y^3[n1  -  y^3{-  +  ot-|n-|  1  ogn-i  +  0(n  )} 

+  ^  +  o(n2))  +  ...] 

W  %  [n-j  -  y-j  ^  n-j  +  0 (ti -^ ) } 

+  y^3{-j2f  n3  +  0(n2)}  +  ...] 


(3.8) 


In  ref.  26  there  is  assumed  to  be  an  inviscid  region  between  the  wall 
layer  discussed  above  and  the  main  part  of  the  streamwise  boundary  layer  which 
has  been  forced  away  from  the  wall  by  the  colliding  crossflow  boundary  lay¬ 
ers.  There  is  an  analogous,  though  not  identical,  solution  to  equations  (3.2) 
in  this  situation.  We  neglect  the  viscous  terms  and  note  that  pfi/p  is 
no  longer  a  constant  in  this  region  but  Is  given  by  (3.5)  though  we  take  the 
same  value  for  K  as  defined  above.  If  V  is  eliminated  between  the  inviscid 
forms  of  (3.2a),  (3.2b)  and  3/2  k(3/2  k  +  l)pg/(pK)  is  denoted  by 
-R(U),  then  they  may  be  solved  for  W  in  terms  of  U  to  give 


W  -  1/2  U 


[R(U)  -  1/4  IT] 


when  R(U)  >  1/4 


[X  -  X  +  MU)][R(U)  -  1/4  U2]1/2  ,3  q) 

-  =  cot  {  —£ - 1 - 0 - }  l3*9J 

which  holds  near  the  wall  where  U  «  0  and  R(U)  *  1,  but 


[X$  -  X  +  fj(U)][l/4  U2  -  R(U)]1/2 
D - 


W  -  1/2  U 


=  coth  { 


)  (3.10) 


when  R (U)  =  1/4U2  which  holds  as  Y  +  ®  where  1/4U2  -  R ( U)  =  1/4K  (3k  +  l)2. 


Here  f-|(U)  is  an  arbitrary  function  of  U,  and  U  is  assumed  to  vary 
monotonically  from  its  value  0  at  the  wall  to  K  '  as  Y  +  ®.  The  solutions 

(3.9),  (3.10)  continue  smoothly  through  the  zero  of  R(U)  -  1/4  U2.  At  Y  *  ®, 

(W  -  1/2  U) (1/4  U2  -  R(U))-1^2  =  -(3k  +  l)/|3k  +  1|  so  is  negative  if  k  > 

-1/3  and  if  the  inviscid  solution  is  to  satisfy  the  boundary  condition  then 
the  right-hand  side  of  (3.10)  must  tend  to  -1  as  U  tends  to  K"^2.  This 

means  that  f^(U)  -*•  as  U  -*•  K"^2  but,  as  shown  below,  f^(U)  >0  as  U  ■* 

0  so  fj(U)  must  vanish  for  some  U  >  0.  As  the  wall  layer  is  approached  W  is 
large  and  positive  while  U  is  proportional  to  X$  -  X  and  is  thus  negligible. 
Thus,  if  k  >  1/3,  W  -  1/2  U  has  to  change  sign  and  this  will  occur  in  the  region 
where  (3.9)  holds  as  the  hyperbolic  cotangent  cannot  vanish.  It  is  not  clear 
in  which  region  f-j(U)  will  vanish  for  nonzero  U.  This  zero  of  f^(U)  in 
addition  to  that  which  occurs  when  U  =  0  will  imply  that  W  is  again  very  large 
and  could  lead  to  an  additional  singularity;  this  time  in  the  interior  of  the 
fluid,  but  further  numerical  work  is  required  to  verify  this.  Our  prototype 
solution  is  at  k  *  -1/2  which  is  in  the  range  -2/3  <  k  <  -1/3  where  the  dif¬ 
ficulty  does  not  occur. 

From  (3.8)  it  can  be  seen  that  as  soon  as  |f-|(U)|  »  (X$  -  X),  W  is  a 
function  of  U  alone,  and  a  solution  of  (3.2)  with  the  viscous  terms  retained 
may  be  found  in  the  form 

u  *  W  +  (Xs  *  +  0{(Xs  '  X)2} 

V  -  Q ' ( X ) U  *  VQ( Y-| )  +  0(XS  -  X)  (3.11) 

W  *  WQ( Y-j )  +  (Xs  -  XJW^Y^  +  0((X$  -  X)2} 
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where  Y1  =  Y  -  Q(X).  Here  U^Y-j),  WQ( Y1 )  are  functions  of  Y^ 
satisfying  UQ  -►  K“1/2,  WQ  -*•  -  3/2  kK"1/2  as  Y]  -*•  «,  and  UQ  ♦  0,  W  ♦  • 
as  Y1  -*•  and  VQ,  V-j ,  may  be  determined  successively.  This  was  also  the 

situation  in  ref.  26  where,  as  here,  the  function  Q(X),  which  is  below  shown 
to  be  large  as  X  -*■  X^,  gives  the  position  of  the  main  part  of  the 
streamwise  boundary  layer. 

The  matching  between  (3.8)  and  (3.9)  follows  closely  that  in  ref.  26. 
Since  (3.6)  is  fully  determined  apart  from  e-|  in  (3.7)  and  an  arbitrary  con¬ 
stant  y.|  at  the  n  =  1  stage,  both  of  which  are  expected  to  be  determined 
by  conditions  for  X  <  X$,  the  matching  serves  merely  to  determine  f^U) 
for  small  U  and  to  verify  the  consistency  of  the  procedure.  When  U  is  small 
but  W  large,  as  occurs  at  the  outer  edge  of  the  wall  layer  as  was  shown  in 
ref.  26  by  an  examination  of  (3.6)  for  n-j  »  1.  it  follows  from  (3.9) 
that 

W%  U/[X$  -  X  +  f^U )]  (3.12) 

and  matching  with  (3.6)  as  in  ref.  26  leads  to 

f(U)  *  -3-^  [1  -  a,(la?  +  1)  0  (z)]  (3.13) 

2U  J/  ice.  y  u 


where  U  ■  4/3  log  8^/U.  In  ref.  26  it  was  possible  to  define  a  meaningful 
stream  function  for  the  inviscid  flow  that  held  for  all  Y.  Here,  however. 


this  is  not  so  but  if  we  wish  we  may  define  \|>  such  that 

.  U 
FY  '  W 


(3.14) 
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when  U  «  1,  W  »  1,  l.e.  when  »  1.  It  also  follows  from  the  form  of  U 
in  (3.8)  that  in  this  region  U  itself  is  a  function  of  \J>  with 


U 


[1  ♦  (- 


3a, 


1  a-j  log  ¥ 

4)  $ 


+  0  (i)]. 


(3.16) 


where  4*  *  4/3[a-|  1  og (  -j/^)],  and  that  f-j(U)  =  f (\|>)  where 


f(*)  = 


A 


24' 


'm 


3a2  ^  a1  log4' 


[i  -  Hr  -  i) 


+  0  (i)] 


(3.17) 


The  advantage  of  the  formulation  in  terms  of  \|>  is  that  Q( X)  may  be  cal¬ 
culated  as  was  the  corresponding  quantity  in  ref.  26.  On  the  assumption  that 
the  inviscid  flow  contains  no  singularities  other  than  as  U  0,  i.e.  there 
is  no  value  of  Y  at  which  f^  (U)  in  (3.10)  changes  sign  when  U  =  0(1),  then 
the  principal  contribution  to  Q(X)  comes  from  the  inner  part  of  the  inviscid 
layer  and 

fi  7/1  28a.  log  y.  _a/7 

Q(X)  *  [I  -  ~g - 373-  +  0  (yj  '  )J  (3.18) 

1  yl 


Thus  the  match  appears  to  be  consistent.  If  K  <  0,  i.e.,  k  >  0  or  k  <  -2/3, 
the  procedure  cannot  be  carried  through  as  yj  then  turns  out  to  be 
negative  and  it  is  not  possible  to  regard  the  situation  as  the  termination  of 
a  boundary  layer  that  exists  for  X  <  Xg.  If  0  >  k  >  -1/3  it  may  be 
that  an  additional  viscous  layer  is  required  in  what  we  have  termed  the 
inviscid  region.  However,  even  in  that  situation  the  match  between  the  wall 
layer  and  the  inviscid  region  is  satisfactory  and  the  additional  layer  will 
affect  only  Q(X)  in  (3.18)  and  the  main  part  of  the  streamwise  boundary  layer 
that  is  now  centered  on  Y  s  Q(X). 

In  the  following  section  we  shall  make  comparisons  between  the  theory  and 
a  numerical  integration  at  k  *  -1/2,  but  before  doing  so  it  is  worth  consid¬ 
ering  briefly  the  effect  of  suction  on  the  present  asymptotic  theory  since 
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suction  is  used  to  initiate  this  computation.  The  change  in  the  boundary  con¬ 
ditions  required  is  that  now  in  (3.3)  V  =  Vw  where  Vyis  a  given  negative 
constant.  The  leading  term  of  the  expansion  (3.6)  is  formally  unaltered  but 
the  compatibility  condition  of  ref.  26  is  now  dominated  by  instead  of  by 
K  so  that  instead  of  (3.7)  we  have 

V  log(y  1/3+  a^)  =  log  (3.19) 

where  a]  =  -21/^1,  a2  is  a  constant  determined  solely  by  the  singularity  and 
$1  a  constant  which  depends  on  the  history  of  the  solution.  The  expansion  (3.6) 

now  proceeds  in  descending  powers  of  yf^73:  thus  whether  there  is  suction  or 
not  the  appropriate  powers  of  log  [B-|/(XS-X)]  are  negative  integers.  If  Vw  >  0 
the  theory  fails.  It  is  possible  that  in  this  case  an  analysis  similar  to 

07 

that  carried  out  by  Catherall,  Stewartson  and  Williams  for  the  separation 
of  a  boundary  layer  with  uniform  external  velocity  and  uniform  injection  will 
be  relevant. 


3.4  Comparison  with  the  Numerical  Solution  When  k  *  -1/2 

The  numerical  solution  when  k  *  -1/2,  some  of  whose  principal  properties 
are  displayed  in  Fig.  3,  may  be  used  to  illustrate  our  theory.  We  make  the 
comparisons  on  the  scaled  quantities  In  (3.12)  which  are  easily  obtained  from 
the  numerical  Integration  of  (2.19)  to  (2.21).  We  define  the  longitudinal  and 
azimuthal  skin  frictions  by 

TX  ■  57  1,0  *"d  -  It  )y,0  <3-20> 

respectively,  and  it  then  follows  as  in  ref.  26  that 
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(3.21) 


^(X)  *  (tJ/3  +a1a2)3a2/4  exp  (-^t473) 
should  be  close  to  (X$  -  X)/^  as  X*  Xg,  and  that 


TX 


9,(X)  =  —  (1  +  0.266t"8/3) 


(3.22) 


should  be  close  to  X$  -  X.  We  may  also  compare  Wmax  which  is  given  by 


the  formula  of  ref.  26  as 


<t> 


(3.23) 


on  identifying  with  t^,  which  is  consistent  with  the  order  of  accuracy  of 
(3.23). 

In  Table  1  we  present  t^,  tx,  ^(X),  Q^X)  and  in  Table  2  present  Wmax 

computed  from  the  differential  equations  and  also  from  (3.23).  For 

convenience  we  retain  the  abscissa  as  x  and  note  the  relationship 

(xc  -  x)/xc  4  (X.  -  X)  when  X  -  X  is  small.  From  ti(X)  we  would  infer  that 

x$  =  0.3786  and  =*  53.8  while  from  ©-j(X)  that  x$  =  0.3793.  The  agreement 

between  the  two  calculations  of  W_!1V  is  also  satisfactory  in  view  of  the 

max 

order  of  magnitude  of  the  error  in  (3.23). 

One  comparison,  possible  in  ref.  26,  but  not  here  because  of  the  too  rapid 
rise  in  the  boundary-layer  thickness  as  x  ♦  x”.  Is  between  the  asymptotic 
estimate  for  the  displacement  thickness  and  the  value  computed  from  the 
differential  equation.  The  significant  function  is 

1/2  11m  (V  +  1  (k  +  1 )YK”1/2> 

Y-k»  c 


Q'(X)  »  K 


(3.24) 


Table  1.  The  Values  of  ta,  tx  from  the  Numerical  Integration 
and  i;i(X),  9i(X)  from  (3.21),  (3.22) 


X 

x4> 

TX 

SlU) 

e  i  ( 

0.20 

2.202 

0.826 

0.01264 

0.387 

0.25 

2.349 

0.669 

0.00853 

0.293 

0.30 

2.564 

0.478 

0.00473 

0.191 

0.35 

2.978 

0.220 

0.00148 

0.074! 

0.3562 

3.074 

0.180 

0.00112 

0.059! 

0.3625 

3.191 

0.137 

0.000796 

0.043( 

0.3656 

3.272 

0.116 

0.000629 

0.035; 

0.3672 

3.312 

0.104 

0.000559 

0.031( 

0.3687 

3.366 

0.0921 

0.000476 

0.027  * 

0.3703 

3.419 

0.0797 

0.000406 

0.023( 

0.3719 

3.492 

0.0676 

0.000328 

0.019! 

Table  2. 

The  Values  of  Wmax  from  the  Numerical  Integration  and 
from  (3.23),  and  of  Q* ( X)  from  (3.25) 

X 

Hnax 

Eq. (3.23) 

Q '  (x) 

0.20 

2.082 

1.784 

9.4' 

0.25 

2.345 

2.463 

13.2! 

0.30 

2.842 

3.130 

22.2! 

0.35 

4.264 

4.905 

66.0' 

0.3562 

4.876 

5.915 

86.3J 

0.3625 

5.327 

6.081 

122.11 

0.3656 

5.774 

6.570 

153.11 

0.3672 

6.050 

6.823 

175.2' 

0.3687 

6.376 

7.177 

204.01 

0.3703 

6.749 

7.534 

245.5' 

0.3719 

7.246 

8.036 

300.3! 

and  according  to  our  theory 

2t^  4a,  log  t  _.n 

Q'(X)  *  — ^  (1  +  o(t“  /3)}  as  (3 

X  t4> 

The  values  of  Q'(X)  from  this  formula  are  also  displayed  in  Table  2. 


3.5  Computation  of  General  Boundary  Layers 

As  with  our  studies  on  leeward  line  of  symmetry,  a  nonsimilar  solution  was 
produced  near  the  vertex  of  the  cone  by  assuming  it  to  be  permeable  and 


sucking  fluid  away  from  the  boundary  layer.  Specifically  we  set 


P  =  < 


+0.5  exp[-(50x)J] 


0  <  x  <  0.05 


x  >  0.05 


(3.26) 


at  n  =  0  for  0  <_  <J>  <_  n  where  4>  sin  6C  =  9  and  &c  =  10°.  Now  a  solution  of 
(2.19)  -  (2.21)  is  possible  at  x  *  0  by  neglecting  the  derivatives  with 
respect  to  x  and  solving  the  two-dimensional  bound ary- layer  equations  which 
result,  using  the  standard  box  method  and  Integrating  from  the  windward  line 
of  symmetry  <|>  *  0.  At  each  new  station  of  x  the  integration  is  carried  out 
in  the  same  way  with  the  suction  velocity  defined  in  (3.26)  until  at  x  *  0.05 
it  is  switched  off  being  then  already  extremely  small.  Further  downstream  the 
standard  box  method  was  modified  as  necessary  by  introducing  the  zig-zag  box 
but  it  is  not  clear  that  much  was  achieved  thereby.  For  the  limit  of  acces¬ 
sibility  of  integrations  from  the  windward  side  is  virtually  coincident  with 
the  vanishing  of  g"  at  n  =  0  and  It  is  only  when  g'  <  0  that  the  zig-zag 
box  makes  for  a  significant  Improvement.  Calculations  from  the  leeward  line 
of  symmetry  are  desirable  to  complete  the  accessibility  line  and  would  need 
the  more  sophisticated  techniques  of  the  zig-zag  and  characteristic  boxes 
since  g'  changes  sign  across  the  boundary  layer  near  4>  ■  n.  However,  certain 
difficulties  emerged  from  the  studies  near  4»  s  tt  which  need  to  be  overcome 
before  this  can  be  Initiated.  We  shall  refer  to  them  below  and  hope  to 
address  them  at  some  future  date. 
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The  first  calculation  carried  out  was  for  k  =  -0.5  where  k  is  defined  by 
(2.23).  This  value  was  chosen  as  being  in  the  range  when  the  leeside  solution 
develops  a  singularity  at  a  finite  value  of  x  but  the  semi-similar  solution, 
valid  in  the  limit  x  -*■  »,  exists  for  all  <j>  <  tt  and  implies  a  boundary- 
layer  collision  at  $  =  it.  The  step  length  in  x  in  the  complete  calcula¬ 
tion  varied  from  0.005  in  the  suction  regime  to  0.025  when  the  solution  was 

• 

almost  semi  similar  (x  <=),  inn  was  0.1  and  in  4>  varied  from  10°  near 
the  windward  line  of  symmetry  to  0.2°  right  near  the  leeward  line.  The 
results  were  straightforward  as  expected  except  near  the  leeward  line  and 
deserve  no  further  comment.  In  this  neighborhood,  however,  an  interesting 
nonuniformity  developed  which  is  illustrated  in  Fig.  4  where  the  computed 

II 

values  of  f  are  displayed  as  a  function  of  <J>,x.  Three  estimates  of 
f"  at  <j>  =  180°  are  shown,  one  obtained  by  integrating  with  A<t>  =  5°,  the 
second  by  integrating  from  <t>  =  175°  with  A<j>  *  1®  and  the  third  by  inte¬ 
grating  from  <j>  =  179®  with  A4>  *  0.2®.  All  three  differ  substantially  and, 
for  x  <  0.375,  from  the  leeside  line  of  symmetry  solution  although  it  could 
be  argued  that  they  are  all  consistent.  The  most  likely  explanation  is  that 
there  is  an  eigenvalue  attached  to  the  leeside  line  of  symmetry  solution  so 
that  when  180-4>  is  small 

f  =  f£(n,x)  +  An(180-<J))n  fAn(n.x)  +  ...  (3.27) 

where  n  is  not  an  integer  and  An  is  arbitrary.  We  have  not  been  able  to 
pursue  this  possibility  in  detail  but  a  preliminary  inspection  suggests  that 
n  j;  0.2.  Some  additional  evidence  In  favor  of  a  nonintegral  eigenvalue  Is 

I  O 

Indirectly  provided  by  Murdock's  studies  on  llne-of -symmetry  solutions  but 
we  should  emphasize  that  none  has  been  found  to  date.  The  consequences  for  an 
Integration  procedure  starting  from  the  leeside  may  be  serious  since  An 
depends  in  some  way  on  the  rest  of  the  solution  and  cannot  be  specified  a 
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Fig.  4 


.  Effect  of  varying  step  size  in<:the  circumferential  direction  on  the 
wall  shear  parameter  of  the  leeward  line  of  symmetry,  k  e  -0.50. 

The  symbols  o — o—  -denote  leeward  solutions  obtained  from  the 

line- of- symmetry  equations.  The  solid  lines  ( - )  denote  solutions 

obtained  from  the  full  equations  with  a#  =  5°,  the  dashed  lines  those 
with  a<|>  *  1°  for  <f>  >  1750,  and  the  line  denoted  by  A  is  obtained  by 
taking  A$  *  0.2°  for  4  >  179°. 
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priori.  It  is  possible  that  and  n  can  be  fixed  when  x  is  small  but  this 
remains  to  be  investigated. 

The  implication  of  our  analytic  studies  is  that  the  cross-flow  velocity  is 
zero  on  the  leeside  line  of  symmetry  for  x  <  xs  (=0.3793)  and  thereafter 
there  is  a  collision  phenomenon  defined  by  g  having  a  nonzero  limit  as  <J>  -► 
it.  It  follows  that  the  cross-velocity,  in  the  two  boundary  layers  which 
spread  around  each  side  of  the  cone  from  the  windward  line  of  symmetry,  are 
nonzero  and  equal  in  magnitude  but  opposite  in  sign  when  they  reunite  at  the 
leeside.  The  consequence  is  that  a  jet  of  fluid  is  emitted  from  the  boundary 
layers  here  although  the  present  theory  does  not  enable  us  to  make  any  firm 
deductions  about  its  subsequent  fate. 

Our  numerical  studies  lend  some  support  to  the  predictions  of  the 

•I 

analysis.  For  example,  when  x  <  xs  the  limit  value  of  g  as  4>  tt"  is  less 
than  10~5  in  magnitude  (typically  it  is  10"1  for  <|>  <  tt)  and  oscillates  in  sign. 
Within  the  limits  of  numerical  error  we  may  infer  that  it  is  zero  at  4  =  ?r.  If 

II 

x  <  x$  the  limit  value  of  gw  soon  becomes  definitely  positive,  of  order 
-4 

10  ,  and  seams  to  be  approaching  the  value  required  by  Fig.  2b  as  x  -►  ®. 

The  next  calculation  was  carried  out  for  k  *  -0.75  when  the  leeside  line- 
of-symmetry  equations  again  have  dual  solutions  and  the  semi-similar  equations 
valid  as  x  +  «,  have  solutions  for  all  |<t»|  <  tt  a  collision  occurring  at 
4»  =  tt.  The  principal  features  of  interest  are  again  in  the  neighborhood 
of  <j>  =  tt.  In  Fig.  5  we  display  the  variation  of  f"  with  x  for  various 
<|>  and  we  see  that  for  4>  <  170°  the  structure  of  f"  Is  entirely  rea- 

™  n 

sonable,  the  limit  solution  being  rapidly  obtained  as  x  Increases  once  the 
suction  Is  removed.  For  4>  =  175°,  the  smooth  trend  downwards  of  f"  is 
disturbed  and  at  the  next  station  of  namely  the  leeside  $  *  180°,  f 

wt 

Increases  and  shows  little  sign  of  becoming  constant  as  x  Increases. 
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Moreover,  there  is  a  semi -similarity  solution  along  <t>  =  180°  which  is  quite 

different  from  the  limit  form  as  <t>  180°.  It  is  possible  that  closer 

agreement  between  the  two  calculations  will  be  obtained  if  A<j>  is  decreased 

substantially  near  <J>  =  180°,  and  also  that  the  leeside  solution  has  an 

associated  eigensolution  at  k  =  -0.75  as  well  as  at  k  =  -0.50. 

The  third  calculation  was  carried  out  for  k  =  -0.82  when  the  semi-similar 

solution  valid  as  x  ■+  ®  separates  at  <t>  »  170°.  Cooke  examined  such  solu- 
21 

tions  for  k  <_  -1.33  and  concluded  that  they  develop  a  singularity  as  the 
separation  point  *  4>s  is  reached,  but  we  see  no  evidence  of  it  when 

II  II 

k  =  -0.82  from  the  behavior  of  f  and  g  .  Even  when  the  step  length  in  (j>  is 

w  w 

reduced  to  1°,  behavior  of  g"  shows  no  sign  of  the  square-root  behavior 

W  ** 

9“w  *  (4>s  "  <t>)1/2  (3.28) 

?8 

predicted  by  Brown's  theory  .  The  integration  of  the  general  equations 
showed  that  even  with  the  use  of  the  zig-zag  box,  great  care  must  be  used  near 
the  separation  line.  Thus,  after  the  suction  is  switched  off,  we  uniformly 
took  Ax  =  0.01  but  either  chose  A4>  *  5°  or  reduced  it  to  A4>  =  1°  for  <)>  >  150°. 
In  the  first  instance,  the  solution  could  be  obtained  right  up  to  4>  *  180° 
with  largely  believable  data  for  <!>  >  165°  even  though  g"w  <  0.  Moreover,  eo, 
the  effective  blowing  velocity  of  the  boundary  layer,  defined  by 

00  CO  CO 

w 

(l-rldn  +  x^  1(1-  f')dn  +  — 2  /  (g1  -  we)dn 
0  0  2ue  0 

CO 

•  irk  1  (s'  ■  "e,dn  (3-29) 

e  o  ' 

shows  little  sign  of  irregular  behavior  in  the  separated  region. 

Once  the  step  length  A4>  is  reduced  to  1°  a  dramatic  change  takes  place 
in  the  solution:  g"w  changes  sign  at  x  *  0.08,  <t>  *  168°,  for  x  =  0.09 
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the  solution  breaks  down  at  <j>  =  168°,  and  by  x  =  0.120  the  separation  line 
has  been  settled  down  to  4>  *  166°.  As  with  the  semi -simi lari ty  solution, 
when  x  =  «,  g"  shows  little  sign  of  any  singular  behavior  but,  on  the 
other  hand,  era  peaks  up  to  137  at  <j>  =  167°,  x  =  0.100.  This  last  result 
is  in  marked  contrast  to  the  behavior  for  <{>  <_  165°  where  it  is  less  than  4. 

The  main  conclusion  from  this  integration  is  that  unless  care  is  taken, 
the  boundary-layer  equations  will  generate  solutions  where  none  should  exist 
and  moreover  the  effect  of  the  boundary  layer  on  the  external  flow,  in  the 
neighborhood  of  separation,  will  be  diminished.  Since  we  know  from  our 
studies  of  two-dimensional  separation  that  it  is  important  to  have  accurate 
estimation  of  the  displacement  thickness  in  order  to  carry  out  effective 
interaction  studies,  this  conclusion  is  a  matter  of  some  concern.  Of  course, 
it  may  be  that  the  effect  is  weaker  in  three-dimensions  than  in  two  but  it 
would  be  premature  to  assume  this  at  the  present  stage.  It  was  hoped  to  use 
the  calculation  here  to  comment  on  the  parabolic  elliptic  method  for  solving 
boundary-layer  problems  on  cones.  In  this  method  diffusion  terms  in  4>  are 
added  to  the  equations  and  the  new  set  solved  by  the  ADI  procedure.  In  the 
event  we  cannot  do  this  at  present  because  we  have  not  obtained  the  correct 
solution  in  any  part  of  the  integration  domain  <J>  >  <j>s .  A  further  study  of 
the  solution  near  the  leeward  line  of  symmetry  seems  called  for  in  light  of 
the  present  work  to  elucidate  the  role  of  the  eigensolution  and  to  mark  out 
the  remainder  of  the  accessibility  boundary  possibly  with  the  use  of  the 
characteristic  box. 

A  final  calculation  was  carried  out  to  test  our  conclusions  about 
separation.  This  was  done  at  k  *  -0.95.  Here  breakdown  first  occurred  in  the 
suction  region  and  by  the  end  of  this  region  separation  was  well  established 

II 

at  $  =  153°.  Again  there  is  little  sign  of  a  singularity  in  g  as  separation 


is  approached  and  while  a  peak  in  ea  was  observed  at  x  =  0.05,  4>  =  153°. 
However,  the  magnitude  of  this  peak  is  only  14  and  much  weaker  than  when 
K  =  -0.82.  Comparative  plots  of  the  variation  of  eo  with  <j>  for 
different  values  of  x  are  shown  in  Fig.  6. 


3.6  Summary 

We  have  shown  that  there  is  strong  evidence  that  a  singularity  can  and 
does  develop  in  the  integration  of  the  boundary-layer  equations  on  the  leeward 
line  of  symmetry  of  a  cone  at  incidence.  Further  we  believe  that  there  is 
convincing  evidence  that  when  it  does  occur  it  heralds  the  start  of  a  collis¬ 
ion  pheomenon  in  which  the  boundary  layers  coming  from  the  windward  line  of 
symmetry  on  the  two  sides  of  the  cone  have  a  nonzero  cross-flow  velocity  dis¬ 
tribution  when  they  meet  at  the  leeside  line  of  symmetry.  A  full  resolution 
of  this  phenomenon  is  not  yet  available  and  almost  certainly  will  require  con¬ 
sideration  outside  the  scope  of  a  hierarchical  boundary-layer  theory. 

The  formal  conditions  for  the  occurrence  of  the  singularity  were  deduced 
by  adapting  the  theory  in  ref.  26  for  the  singularity  in  the  boundary  layer  on 
the  inside  line  of  symmetry  of  a  curved  pipe  and  are  that  0  >  k  >  -2/3  what¬ 
ever  be  the  Mach  number  and  prescribed  wall  temperature.  However,  it  is  known 
that  similarity  solutions  of  the  governing  equations  exist  when  k  >  k^ 
and  k >  k  >  -1  where  k1*(<0)  and  kg*  depend  on  Mg  and  the  wall  temperature 
although  k-|*  >  kg*  >  -2/3.  Our  numerical  solutions  indicate  that  if  the  initial 
profile  is  arbitrary,  these  similarity  solutions  will  appear  as  limits  of  more 
general  solutions  when  x  -*■ 00  and  that  only  if  kg*  <  k  <  k-j*  does  the 
singularity  occur. 


Fig.  6.  Variation  of  e„,  the  reduced  normal  velocity  at  the  outer  edge  of  the 

boundary  layer,  with  4  *  e/slns.  for  k  =  -0.95  and  x  *  0.025,  0.05,  0.1. 
It  Is  noted  that  similarity  Is  virtually  achieved  as  soon  as  the  suction 
Is  removed. 
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